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, Abstract. We study the dynamics of a single electron spin in a double quantum dot 

(DQD) and its readout via a quantum point contact (QPC). We model the system 
microscopically and derive rate equations for the reduced electron density matrix of 
the DQD. Two cases with one and two electrons in the DQD are studied. In the 
one-electron case, with different Zceman splittings in the two dots, the electron spin 
i-^j , states are distinctly characterized by a constant and an oscillatory current through 

£^ ' the QPC. In the two-electron case, the readout of the spin state of the electron in 

one of the dots called the qubit dot is essentially similar after considering hyperfine 
interactions between the electrons and the nuclear spins of the host materials and a 
uniform magnetic field applied to the DQD. Moreover, to ensure that an electron is 
properly injected into the qubit dot, we propose to determine the success of the electron 
£■ — ' injection from the variations of the QPC current after applying an oscillating magnetic 

■ field to the qubit dot. 
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A single or a small number of electron spins confined in a semiconductor quantum 
dot (QD) has become a subject of considerable interest, partly motivated by potential 
applications in quantum information processing. Because a electron spin in a QD can 
have a relatively long decoherence time, it is a promising candidate for realizing a qubit 
[U [2] , the basic unit of a quantum computer. For a single electron spin, the dephasing 
time is found in theory to be /is in both GaAs [31 H] and InAs dots [5]. Experimental 
results show that the ensemble dephasing time of an electron spin has the order of ns 
P, [7]. Moreover, the spin relaxation time in a large GaAs dot is found to be about 
1 — 10 /jLS at a moderately low temperature of 10K [8]. Indeed, recent experiments have 
demonstrated that spins in QDs can be used to carry quantum information [9l [TOl HT]. 
For both applications in quantum computing and fundamental research, the readout of 
qubit states based on electron spin is a centrally important issue [12]. However, due to 
the weak magnetic moment associated with the electron spin, it is difficult to directly 
measure the electron spin states. A possible solution is to correlate the spin states to 
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charge states, and the measurement of the charge on the dot will provide information 
about the original spin states j2]. This can be implemented using a quantum point 
contact (QPC), which is a charge detector and can be used to determine the number 
variation of the electrons confined in the QD. 

Recently, the readout of electron spin states in a QD has been realized using such 
a spin-charge conversion [TOj [13]. For the experiment in Ref. [10], a QD is connected 
to an electron reservoir. Applying an external magnetic field, gate voltages are applied 
so that the electron confined in the dot can tunnel to the reservoir if its spin is down. 
(A spin-up electron cannot tunnel in this case). A nearby QPC is used to detect the 
electron number variation in this QD and can determine the electron spin state in the 
QD. Also, Engel et al. [H] proposed various implementations of the readout process 
based on a double quantum dot (DQD). Barrett and Stace [15] proposed a electron 
spin readout approach using a microwave field and an inhomogeneous Zeeman splitting 
across the DQD. 

In the present paper, we study two implementations for reading-out electron spin 
states based on a DQD coupled to a QPC. Also, we explain the effects of static and 
oscillating magnetic fields on the electron spin states. The first implementation involves 
a single electron in the DQD. The readout of the spin states is based on the difference 
Zeeman splittings in the two QDs. In the second implementation, two electrons are 
allowed in the DQD. The Pauli exclusion principle and hyperfine interactions between 
the electrons in the DQD and the nuclear spins in the host materials enable the 
readout of the electron spin states. These are interesting examples for implementing 
readout of the electron spin states. A potential advantage of our proposal is that the 
readout manipulation can easily be switched on (off) by decreasing (increasing) the 
tunneling barrier between the two dots through varying the gate voltages. Thus, the 
readout process can be implemented only when needed. This is important in quantum 
information processing. To understand the underlying physics from a microscopic point 
of view, we derive a set of rate equations describing the electron dynamics of the DQD 
system. Based on these rate equations, we calculate the QPC current and illustrate that 
the QPC current behaves differently for different spin states. 

The paper is organized as follows. In Sec. 2, we model the system when only one 
electron is confined in the DQD. A set of Bloch-type rate equations are derived to 
describe the detailed measurement processes for the electron spin states in the qubit 
dot. In Sec. 3, we study the measurement of the electron spin states in the qubit dot 
when two electrons are confined in the DQD. Sec. 4 is the conclusion. 

2. Readout of single electron spin: one electron in DQD 

2.1. Theoretical model 

We first discuss a scheme to detect the electron spin states in the case with only one 
electron confined in the DQD. As schematically shown in Figure. [H the whole system 
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Figure 1. (Color online). Schematic diagram of a double quantum dot (DQD) and 
a quantum point contact (QPC) with only one electron in the DQD. The left (qubit) 
dot of the DQD is coupled to a right (reference) dot via hopping. The nearby QPC 
is used as a detector measuring the number variation of electron in the reference dot. 
An energy-level detuning of the electron spin states is generated using two external 
magnetic fields in the two dots, (a) Gate voltages are adjusted to keep Eli = E^, 
so that the hopping of the spin-up electron between the left and right dot is allowed. 
Moreover, the energy-level detuning for the spin-down electron is much larger than 
the hoping strength, i.e., Er± — Eli ^o> ano - nence the hopping of the spin-down 
electron between the two dots is forbidden, (b) The hopping blockade for the spin- 
down electron is lifted by applying a transverse magnetic field B(t) = B£ cos(cj c t) to 
the left dot, which flips the electron spin. 

consists of a DQD and a QPC. The left dot is used as a qubit dot, in which the electron 
spin is expected to be readout. The right dot is used as a reference dot. The QPC 
is capacitively coupled to the right dot and serves as a readout device. The electron 
number variation in the right dot induces a change in the barrier in the QPC. This 
leads to a variation of the current through the QPC, which can be used to indicate the 
occupation of the right dot [T6] . 

The Hamiltonian of the whole system is given by 

H = H D q D + Hqpq + Hmt + H rl , (1) 

with 

#DQD=X) E i* c l c ic + tt {c j La C Ra + 4 CT C L(J ), 
in o 

(2) 

where i = L, R denote the left and right dots, and c\ a (c^) is the creation (annihilation) 
operator of electron with spin o in the z'th QD. Q denotes the hopping amplitude 
between the two dots and here it is assumed to be spin- independent. We have denoted 
the energy levels in the z-th dot by -E^Q) = E % ~F \^4-> with = g[i B B*, where Ei is 
the orbital energy level of the QD and is an externally applied magnetic field in the 
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i-th dot along the z direction. Here, g is the effective gyromagnetic factor and \ib is the 
Bohr magneton. We have chosen the unit h — 1. The Hamiltonian of the QPC reads 

-f^QPC = ^akO'ak a ak + [Ak a rk + a tk a lk) > (3) 

afc irfc 

where a^ fc (a a fc) is the creation (annihilation) operator of an electron with momentum k 
in reservoir a (a = l,r). H mt gives the electrostatic interaction between the DQD and 
the QPC: 

H in t = ^irCR a c Ra (aj k a rk + a\ k ai k ). (4) 

Irkcr 

An electron spin resonance (ESR) magnetic field is applied in the x direction at the left 
dot, leading to a term 

H Ti = A x .(t)(c+c Li + H.c.) (5) 

with A x (t) = ^gfisB^ cos(u c t). This ESR magnetic field generates spin flipping when 
it is resonant with the Zeeman splitting on the left dot, i.e., u c = g[iBB z L . 

The spin-up and spin-down states | \t) and | [l) in the left dot constitute the basis 
states of a qubit. The right dot works as a reference dot, the electron occupation of 
which is measured by the nearby QPC. Energy detuning for the spin-up (down) electron 
is e-fm = -Ertq) — -Extq) . We assume that the Zeeman splittings A z are different in the 
two dots. This can be realized, e.g., by applying a micro-size permanent magnet near 
one dot of the DQD [T7j. This leads to a difference e-j- — = A Z R — A Z L in the energy level 
splittings for the spin-up and spin-down electrons. In our consideration, gate voltages 
are adjusted to keep £j ~ 0, so that a spin-up electron can hop back and forth between 
the two dots. Furthermore, we also assume ^> e^Qq, so that hopping is forbidden 
for spin-down electron. However, this spin blockade can be lifted by an ESR magnetic 
field. Here, in the one-electron case, the effects of the nuclear magnetic fields in the two 
dots are neglected. This is because the Zeeman splitting in each dot is much larger than 
the nuclear field in the x and y directions. Moreover, the z component of the nuclear 
field only shifts the energy level and can be included in the Zeeman splitting. 

The physical picture of the electron spin readout is as follows. We first inject an 
electron with either up or down spin into the qubit (left) dot. An initially spin-up 
electron in the qubit dot can hop into the reference dot. This will lead to a variation 
of the current through the QPC. In contrast, for an initially spin-down electron in 
the qubit dot, it will remain stationary because £j_ 3> fio- As a result, no variation 
in the QPC current occurs. Therefore, one can determine the initial electron spin 
state based on the variation of the current through the QPC. However, in practical 
experiments, the injection of electrons into the DQD may not be always successful. 
Without any electron in the DQD, there is also no variation in the QPC current. Thus, 
this simple implementation cannot distinguish between the cases with zero or one spin- 
down electron. To solve this problem, as will be shown below, one can apply an ESR 
magnetic field in the qubit dot. The ESR magnetic field induces spin flipping in the left 
dot. If there is a spin-down electron, it can be converted to the spin-up state by the ESR 
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field and then hop onto the right dot. Therefore, a current variation will be observed in 
the QPC. In contrast, the QPC current will remain unchanged in the zero-electron case 
even in the presence of the ESR field. 



2.2. Bloch-type rate equation 

To describe the physical processes quantitatively, we derive a set of Bloch-type rate 
equations for the reduced density matrix a(t) of the DQD system. Following Gurvitz 
et al.[T8l [T9] , we write the wave function of the whole system in the occupation 
representation as 

!*(*)> = E [M*)cL + M*)4, 

a 

+ E h Lair{t)c ] Lu a\a l + bRairi^c^alai 

Ir Ir 

+ E bLall'rr'{t)A rT a l a r' a l a l l 
l<l',r<r' 

+ E 6 arfi'rr'(*) c ftj4 a J' a » a «' + " " -]l )> 
l<l',r<r' 

(6) 

where bj(t), j = La, Ra, Lair, Ralr, ... are the time- dependent probability amplitudes to 
find the system in the corresponding states. For example, bLairit) denotes the probability 
amplitude for the state with an electron having tunnelled through the QPC barrier (from 
the left reservoir to the right one) at time t, and an extra electron with spin a staying 
in the left dot. The vacuum state |0) corresponds to the state where there is no extra 
electron in the DQD and all the energy levels up to the Fermi energies fiL and fiR of 
the two reservoirs of the QPC are occupied by electrons. 

The relevant electron states of the DQD span a four-dimensional Hilbert space. 
We adopt the notations |1) = | |l) and |2) = \ l L ) for the left-dot states, as well 
as |3) = | fi?) and |4) = | Ir) for the right-dot states. A diagonal element aa(t) 
(i = 1,2,3,4) of the reduced matrix represents the occupation probability of the state 
while an off-diagonal element aij(t) characterizes the quantum coherence. Each an 
is further given by 

a ii {t) = ag ) + aS ) +a£ ) + ..., (7) 

where a\^ is the probability that the DQD is at state \i) after n electrons have tunnelled 
from the left reservoir of the QPC to the right one. In this notation, we have, for example 



M*)l 2 , «8 = EI^)I 2 , <i j = E \b LW rA\ .... (8) 

Ir l<l',r<r' 



The current flowing through the QPC is 

dN(t) 



Qpc(*) = e^-A (9) 
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where N(t) is the number of electrons transported to the right reservoir of the QPC at 
time t. Accordingly, we have 

W*) = E^i n) W- (io) 

n,i 

Substituting the many-body wave function of the whole system into the Schrodinger 
equation i\ip(t)) = H\i[)(t)), one gets a set of differential equations for the probability 
amplitudes bj(t). In the nonequilibrium transport in the QPC with a large voltage 
bias, following Refs. [19] and [20], the Bloch-type rate equations for the reduced density 
matrix a(t) of the DQD are derived by integrating the degrees of freedom of the QPC 
reservoirs. By summing &( n \i) over n, the rate equations for the diagonal elements are 
given by 

&n(t) = ifi (>i3 - cr 3 i) + iA x (t)(a 12 - cr 21 ), 
o- 22 (t) = A)(>24 - o- 42 ) + iA x (t)(a 2 i - cr 12 ), 
cr 33 (t) = in (a 31 - a 13 ), 

a u (t) = itt (a 42 - a 24 ), (11) 
The rate equations for the off-diagonal elements are 

&i2(t) = i{E L[ - E L1 )cr 12 - ^0^32 + A)f"i4 + iA x (t)(a u - a 22 ), 
o-i3 (t) 

&u(t) 

023 (*) 

024 (*) 

<7u(t) = i(E Ri - E m )a M + ifio(T 32 - i^o^u, 

(12) 

where T d = (VD' — vT>) is the dephasing rate induced by the QPC detector [T9] . 
Here we have defined 

D = 27rp L p R n 2 V d , D' = 2np L p R n' 2 V d , (13) 

and 

^v d {v + n- 2 )' (14) 

with A = 2np L p R Q'Q QVd. In Eq. f fT2"j) , the terms proportional to % are due to the 
inclusion of higher-order terms of O(fl 2 fl /V d 2 ) [20]. Also, we assume that the tunneling 
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Figure 2. (Color online). Time evolution of the electron occupation probability in 
the right dot for (a) spin-up | ti) an d (b) spin-down | electron states in the left 
dot. (c) and (d) Time evolution of the QPC currents corresponding to (a) and (b), 
respectively. We have set the parameters as Slo = 0.25 /xeV, x = 0.0025 ^cV and 
T d = 60 MHz. 



couplings depend weakly on the energy, so that Qi r (Ei, E r ) = Q, and Qi r + SQi r = Q'. 
Vd = hl — hr is the voltage bias applied on the QPC and Pl(pr) is the density of states 
in the left (right) reservoir of the QPC. From Eq. (fl2|) . one can see that characterizes 
the exponential damping of the off-diagonal density matrix elements. Now, the QPC 
current is given by 

/qpc(*) = h[vii(t) + a 22 (t)} + h[(T 33 (t) + a M (t)}, (15) 

where Iq = D is the current flowing through the QPC when the right dot of the DQD is 
empty, while Ii = D' is the QPC current when the right dot is occupied by one electron. 
Since Iq ^ I\ in general, one can determine the electron occupation of the right dot from 
the variation of the QPC current. 

We first consider the case without an ESR oscillating magnetic field, i.e., A x (t) = 0. 
Using Eqs. ( |TT|) and ( fl2l) . one can numerically calculate the occupation probabilities an, 
i = 1 to 4. A typical value for the hopping coupling between the two dots in experiments 
is Qq = 0.25 /xev[21]. We have taken parameters so that the initial current of the QPC is 
Io = 1.5 nA if the right dot is empty [TP] , while it equals l\ = 1 nA if there is an electron 
in the right dot. First, consider the case that a spin-up electron is injected into the left 
dot. Figure [2]^a) shows the calculated occupation probability of the electron in the 
right dot. The corresponding current flowing through the QPC is given in Figure. [2](c). 
It shows that the current /qpc starts from the initial value Iq, and then decreases, 
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Figure 3. (Color online). (a)Time evolution of the occupation probability in the spin- 
up state | ti?) in the right dot in the presence of an oscillating magnetic field. Initially, 
the electron in the left dot is in the spin-down state | (b)The corresponding time 
evolution of the QPC current. We have taken = 0.3 fJ.eV, and tu = 0.5 /LteV. 



oscillates and finally converges to a value other than Iq. In addition, oscillations in both 
the occupation probability and the QPC current are observed. These results from the 
fact that the spin-up electron can tunnel back and forth between the dots. In contrast, 
if the electron injected into the left dot is spin-down, it cannot hop into the right dot 
because e± 3> fi . The electron occupation probability in the right dot is hence zero [see 
Figure. EHb)]. The down spin is also reflected in Figure. E](d), where the QPC current 
remains unchanged. Accordingly, one can distinguish between the two initial electron 
spin states from the variation of the QPC current. In short, if the QPC current decreases 
from its initial value, the initial spin state is spin up. Alternatively, if the initial spin 
state is spin-down in the left dot, the QPC current remains unchanged. 

We have repeated the calculation by considering an additional ESR magnetic field. 
Without such a field, we cannot distinguish between the case with no electron in the 
left dot from that with a spin-down electron as discussed in Sec. 2.1. Both give rise to 
no QPC current variation. In the presence of an ESR field on the left dot (i.e., ^ 0), 
for the zero-electron case, the QPC current remains unchanged. However, for an initial 
spin-down electron in the left dot it can flip to the spin-up state induced 

by the ESR oscillating magnetic field. In contrast, the spin-up electron can hop into 
the right dot [see Figures. E3(a)]. This leads to a variation in the QPC current [see 
Figures. [3](b)]. Therefore, these two cases can now be distinguished. 

3. Readout of single electron spin: two electrons in DQD 

3.1. Theoretical model 

We now study the readout of the electron spin states in the left (qubit) dot assuming 
that an additional electron initially occupies the right (reference) dot (see Figure. H]), as 
in a recent experiment [22]. We further assume that the gate voltages of the dots are 
tuned so that no two electrons can simultaneously stay in the qubit dot. Thus, the two 
relevant occupation configurations correspond to two electrons in the right dot or one 
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Figure 4. (Color online). Schematic diagram of a double quantum dot (DQD) and a 
quantum point contact (QPC) with two electrons in the DQD. One spin-up electron 
is initially kept in reference dot by properly adjusting the gate voltages, (a) A spin- 
down electron in the left dot can always hop into the right dot after taking hyperfine 
interactions into account, (b) Transport of a spin-up electron is forbidden due to Pauli 
exclusion. 



electron in each dot. The total Hamiltonian is 

H = Hqpc + H D q D + H int , (16) 

where Hqpq is the Hamiltonian of the QPC in the two-electron case, which has the same 
form as Eq. ([3j), but with Qi r replaced by fl' lr . The Hamiltonian of the isolated DQD 
system after considering both inter- and intradot Coulomb interactions now becomes 

Hdqd — H + H spin , (17) 

where 

H = J2 Ei4 a C ia + n Y,( C LaCRa + H.C.) 

la a 

+ Uin^nn + U LR n La n Ru , . (18) 

i era' 

In the absence of a net nuclear polarization, randomly oriented and fluctuating 
nuclear spins in the host materials give rise to effective magnetic fields Bnl and Bnr 
in the left and right dot, respectively. They results from different local environments 
for the electrons in the respective dots. However, nuclear fields change with a nuclear 
spin relaxation time scale of the order 1 s, which is much longer than any time scales 
characterizing the transport processes of electron. These nuclear effective fields can thus 
be regarded as static fields in our discussion [231 12H 125] • Therefore, we can describe the 
influence of the magnetic fields on the electron spins in the DQD by 

H S pi n — g^B^NL ■ Sl + g^B^NR ' Sr 

+ 9 fi B BUSl + S Z R ) + g^ B Bl cos^K, 

(19) 
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where Sl and Sr correspond to the electron spin in the left and right dots, respectively. 
The third term in Eq. ffl9|) is the Zeeman splitting caused by an external perpendicular 
field. The last term is an ESR oscillating magnetic field in the x direction. In the 
present discussion, we assume that the ESR oscillating magnetic field is only applied on 
the qubit dot. 

The relevant electronic states for the DQD span a five dimensional Hilbert space. 
The basis set consists of double-dot triplets |1) = \T + ) = | Tz,Tii)> 1 2) = = | IlIr), and 
|3> ee |To) = ^=(|TxIa> + |UW), double-dot singlet |4> = \S d ) = ±(\UIr) -\IlU)), 
and single-dot singlet |5) = \Ss) = Tfliii) ~~ I I i?Ti?) ) • Single-dot triplet states are 
excluded due to their much higher orbital energies [26l [H] . In this representation, -£/dqd 
is rewritten as 

#DQD = ^il Z ')( Z 'l 
1=1,2,3,4,5 

+ 9 -^[(B: + <BJf)|3><l| + (B: - iBJf)|3><2| + H.c] 

+ 9 -^[{~B x d - iBj{)|4><l| + (B* d - iBl)\A){2\ + H.c] 

+ fi (|4)(5| + |5)(4|) + ^^(|3)(4| + |4)(3|) 

+ ^ cos(u; c t)[|3)(l| + |3)(2| - |4)(1| + |4)(2| + H.c], 

(20) 

where B d = |(B NL - B NR ), B s = |(B NL + B NR ) + B z ext % and fti = ^gn B B x L . We have 
also introduced energy levels given by 

#1,2 = ^3 T 9HbB z s , E 3A = E l + E r + U lr , (21) 

and 

E 5 = 2E R + U R . (22) 

A critical step in the readout is the hopping to the right dot, where there is a non-zero 
Coulomb energy barrier 

A = E 5 - E A = U R - U LR - {E L - E R ) , (23) 

for the second electron at the right dot if the intra-dot repulsion U R dominates. 
The interaction Hamiltonian between the DQD and the QPC is 

flint = J2 6 ^'ir cj Ra c R^ c m c m(a]k a ik + al k a rk ), (24) 

Irkcr 

In the singlet-triplet representation, it can be written as 

H mt = Y / Sn'\5)(5\(ala rk + a\ k a lk ). (25) 

Irk 

Similar to the one-electron case, the hopping amplitude Q' lr of the QPC and its change 
5Q' lr by either adding or removing an electron in the right dot are assumed to be energy- 
independent, so that Q' lr (Ei,E r ) = Q' and Q! lr + 5D,' lr = Q". As discussed in the one- 
electron case, the left QD is a qubit dot and the right dot is a reference dot. The nearby 
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QPC works as a detector to measure the number change of electrons in the reference dot. 
When there is one electron staying in the reference dot, the current flowing through the 
detector is I\. For double occupancy in the reference dot, the detector current becomes 
J 2 , where I 2 < h because of the increased QPC barrier induced by the additional 
electron. As a result, the variation of the electron number in the reference dot can be 
reliably detected from the current change of the detector. 

We first briefly discuss the readout processes of the qubit states. We assume that 
the electron that is always kept in the right dot is spin-up. This can be realized by 
injecting an unpolarized electron and wait until a time interval much longer than the 
typical relaxation time of a single electron spin. It will relax to its ground spin-up state 
due to its coupling to the outside environments. An additional electron is then injected 
into the qubit dot and its spin state will be readout. For a spin-down qubit electron, the 
initial total spin in the z direction is S z = 0. The DQD system takes either the double 
electron state |3) or |4) with an equal probability of 1/2. If the state taken is |4), the 
electron can directly hop onto the right dot. This hopping is described by the Qq term 
in Eq. (1201) . Otherwise, if it is |3), the electron can transit to state |4), as allowed by 
the B z d term in Eq. ( 1201) . Then hopping into the right dot becomes possible similarly. 
Therefore, the qubit electron can always hop onto the reference dot on the right, leading 
to the state |5) for S z = 0. Due to this hopping, the QPC current changes from ix to 
another value and this indicates the spin-down qubit state. 

In contrast, for a spin- up qubit electron leading to a total spin component S z — 1, 
it form the triplet state |1) with the spin-up electron in the reference dot. Because of 



the application of a large external magnetic field B* xt 3> y(B%), this triplet state |1) 

is far away in energy from other states. Thus, it is decoupled to states |3) or |4) and 
the electron in the left dot cannot hop to the right dot. Therefore, the current flowing 
through the QPC remains constant at I\ and this indicates that the initial qubit state 
is spin-up. 

3.2. Bloch-type rate equation 

To reveal the quantum dynamics of electron states in the DQD system, we derive a set 
of Bloch-type rate equations for the reduced density matrix a{t) of the DQD, also using 
the technique developed by Gurvitz et al[18]. We assume the high Zeeman splitting 



limit, i.e., B* xt 3> J (B%), in order to suppress the effect of the nuclear fields. Spin flips 
caused by hyperfine interactions are then negligible. The many-body wave function 
\^{t)) of the whole system in the singlet-triplet basis is given by 





*(*)>= E Mt)4 + E bur(t)c+ r alai 



i=l,2,3,4,5 Ir 



+ E onvrr'{t)c[a\a r ,aia v + . . . |0). 



l<l',r<r' 



(26) 
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where |0) is the vacuum state and bj(t) is the time-dependent probability amplitudes of 
the corresponding state For example, when j — ilr, with i — 1,2, .. or 5, bj(t) is the 
probability amplitude of the state with the DQD system at state \i) while one electron 
has already passed through the QPC at time t. In addition, we have used cj (q), which 
denotes the creation (annihilation) operator for state \i) in the DQD system. 

Substituting the wave function [Eq. (126)) ] into the Schrodinger equation 

i\^f(t)) = H\i£(t)), and tracing over the reservoir states of the QPC, we obtain a set 
of Bloch-type rate equations for the reduced density matrix a(t) of the DQD system: 



and 



on (t 

033 (* 

055 

012 
0"13(* 

014 (t 
015 

023^ 
024 (* 

025 

034 (t 

035 (* 

045 (t) 



Z.D_a" 41 — iS-(J 3 \ — Z-D + 0"i4 + fiS' + C7i3, 

- iD + a± 2 - iS + a 32 + iD^a 24 + iS^a 23 , 
ifl d (a u - cr 43 ) + iS + a 32 + iS-O^i - iS-a 23 - iS + a X3 , 
ifl d (a 43 - a 34 ) + iD + a 42 - iD_a A1 + i^o(cr 45 - cr 54 ) 

- iD_a 2 4 + iD + a 14 , 
^0(054 - 045), 

i(E 2 - Ei)<j 12 + iD-a u + iS_a 13 + iD_a 42 - iS_a 32 , 
i(E 3 - £7i)t7i3 + iZL.043 - zSLtTaa + iS+aia 

+ iS'_(Jii + «^(Ti4, 

i(£ 4 - £1)014 + iD-{a u - - iS-au + 

+ ifi O 015 + ^^13, 

i(E 5 - Ex)a X h + i& (Tu + iD-atf - iS^a 35 - 7, T 'd a ^ 

1 , 

- ~ 2 X *u. 

i(E 3 - E 2 )a 23 - iD + a 43 - iS + {a 33 - o 22 ) 
+ iS-a 21 + iVt d a 2A , 

i(E 4 - £2)024 - iD + (au - 022) ~ ^+034 
+ i{l a 25 + itl z d cr 23 , 

i(E 5 

1 , 

- 2* °"24, 

- iS-a-u - iS + a 1A - iD_a 31 + iD + a 32 - iVL d {a 44 - a 33 ) 
i(E 5 - £3)035 + itt (T 3 4 - iS^a 25 - iS + a 15 - i£l d a 45 

V 1 , 



iD_a 21 



E 2 )a 25 + iQ a 24 - iD + a 45 - iS + a 35 - -IV25 



i(E 5 - £4)0-45 + ^0(0-44 
- itt z d (T 35 



a-55) - iD_a 25 + iD + a 15 



(27) 
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Here the detector-induced dephasing rate is T' d = (V D" — s/D'^j , with 

D" = 2np L p R n" 2 V d , D' = 2np L p R n ,2 V d . (21 
Also, we have defined 



d± (t) = + fii (t ) , 5± (t ) = nf + fii (t) , (29) 



where 



A' = 27Tp L p R tt"tt tt'V d , fii(t) = fi a cos(cj c t), 

«J = ^(BS ± iBS), 0? = 9 -^(B: ± LBJf). (30) 
The QPC current is given by 

I(t) = h[a n (t) + a 22 (t) + <7 33 (t) + (744(t)] + /2<755(t). ( 31 ) 

where Ji (J2) is the stationary current through the QPC when the right dot is occupied 
by one electron (two electrons). 

3.3. Results and analysis 

We have numerically integrated the rate equations and obtained the time-dependent 
density matrix elements. As discussed in Sec. 3.1, since there is always a spin-up electron 
in the reference dot, the injection of a spin- up electron into the qubit dot forms a double- 
dot triplet state |1) = | ]l]r) in the DQD system. In contrast, if the injected electron 
is spin-down, the DQD system initially takes the state | [l]r)- Thus, after injecting an 
electron into the left dot, the possible experimental initial states of the DQD system are 
|1) and I |lT-r)- In order to show how the current through the QPC changes for different 
initial states of the DQD system, we assume that the DQD system initially takes the 
state |1) or | Il\r)- 

The initial state | [l\r) is a superposition of the double-dot triplet state |3) and 
the double-dot singlet state |4), i.e., 

|UTh> = ^(|3>-|4)). (32) 

Here the state |4) is coupled to the single-dot singlet state |5) directly via hopping 
coupling, while the state |3) is coupled to |5) via the intermediate state |4) (where the 
transition from |3) to |4) is induced by the B z d term). To reveal the contributions by 
different components, we also take the state |3) or |4) as the initial state to study the 
time evolution of the current through the QPC. 

In our numerical calculations regarding the initial state | |lTr)j we rewrite the rate 
equations (1271) and (1281) in the occupation representation defined by the basis states 
\a), \b), \c), \d), and |e), where 

|a) = |TiTfl> = |l>, |&) = |UU> = |2>, 
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Figure 5. (Color online). Time evolution of the occupation probability 055 (i) for 
the single-dot singlet state \Ss) with four different initial conditions: (a) oii(0) = 1, 
(c) 033(0) = 1, (e) 044(0) = 1 and (g) <Jdd{0) = 1, where the hyperfine interaction 
is not considered, (b), (d), (f) and (h) The corresponding QPC currents. Here we 
have chosen the following parameters: B^' z = (0,0,0) mT, = (0,0,0) mT, 

uj c = 3.75 AieV, A = 0.25 /icV, n = 0.25 /xeV, x' = 0.0025 /ueV, fti = 0.375 fj,eV, and 
T' = 60 MHz. 



|c> = I TlW = ^(|3> + |4», \d) = I UU) = ^(|3> - |4», 
|e) = -L(|t Jl U>-UHTfl» = |5>. 



(33) 



With these new basis states, one can express = (i\cr\j) = 1 to 5) using 

<V = {lA^W) (/-*> v = a,b, c, d, and e), e.g., 

1 , , 1 , 



ai3 = v^I a °" 23 = bc °~ bd 



(34) 



In this way, we can transform equations (12 7p and into the rate equations in the 
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Figure 6. (Color online). Time evolution of the occupation probability a^(t) of the 
single-dot singlet state \Ss) for four different initial conditions: (a) cth(0) = 1, (c) 
033(0) = 1, (e) (744(0) = 1 and (g) <Jdd(0) = 1, where the hyperfine interaction is 
included, (b), (d), (f) and (h) The corresponding QPC currents for these four different 
initial conditions. The nuclear magnetic fields are chosen to be B^' 2 = (—2, 1, 3) mT, 
and Bjjr ~' z = (—1, 2, 0) mT. The other parameters are the same as in Figure. [5] 



occupation representation. 

We first consider the case without hyperfine interactions, i.e., B NL (r) = 0. In this 
limit, the coupling between the states |3) and |4) vanishes, (i) If initially the DQD takes 
the double-dot triplet state |1), the system will not evolve into other states due to the 
large Zeeman splitting. In this case, the electron in the left dot does not hop into the 
right one and the current through the QPC does not change, as shown in Figures. [5]^a) 
and[5](b). (ii) Alternatively, for an initial state |3), the DQD system will also remain 
at this state because |3) does not couple with any other states. Similar to the case of 
the initial state |1), the occupation probability of the single-dot singlet state |5) is zero 
and the current through the QPC also remains unchanged [see Figures. [5]^c) and[5](d)]. 
(iii) In contrast, as shown in Figures. [5](e) and[5](f), if the DQD system initially stays 
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Figure 7. (Color online), (a) Time evolution of the occupation probability of the 
single-dot singlet state in the presence of an ESR magnetic field for the initial state 
|1) (see the text), (b) The corresponding QPC current. 



at 1 4), it couples with the single-dot singlet state due to the hopping coupling between 
the two dots. This gives rise to nonzero occupation probability for the single-dot singlet 
state and a variable current through the QPC. (iv) The results shown in Figures. 0(g) 
and El^h) look like a combination of the results in both (ii) and (iii). This is because 
the initial state \d) = | IlTr) is a superposition of the states |3) and |4) [cf. Eq. (|32l)]. 
Moreover, only the state |4) contributes to the variations of both the probability of the 
state 1 5) and the current through the QPC. 

Moreover, it is shown in Figures \5Hjo) and [5](d) that the two cases with initial 
states |1) and |3) are indistinguishable in measuring the electron spin. This is due to 
neglecting the hyperfine interactions. When they are included, these two cases become 
distinguishable (cf. Figure. [6]). 

For an initial state |1) or |4), the results are similar to those without the hyperfine 
interactions. This can be clearly seen by comparing Figures. [6](a) and [6](b) with 
Figures. [5](a) and Mjo) for initial state |1), and similarly comparing Figures. EJ^e) and 
[6](f ) with Figures. [5](e) and MJ) for initial state |4). In contrast, for initial state |3), 
because |3) and |4) are degenerate, hyperfine interactions are able to provide significant 
couplings. Moreover, state |4) is also coupled to |5) via hopping. Thus, from initial state 
1 3), the system can finally evolve to |5). Indeed, this is reflected in the variations of both 
the occupation probability of state |5) and the QPC current [Comparing Figures. M^c) 
and[6](d) with Figures. M^c) and[5](d)]. For the initial state \d) = | |lT_r), the probability 
of the state |5) and the QPC current [shown in Figures E](g) andEJjh)] also look like a 
combination of the results for both the initial states |3) and |4) [shown in Figures [6](c)- 
E(f)], similar to the case without hyperfine interaction in the DQD. 

The last issue to be addressed is to determine if the electron is successfully injected 
into the left dot. If the injection fails, the QPC current is always I±. The result is the 
same as that with a spin-up electron injected into the qubit dot. To distinguish between 
these two cases, we can apply a transverse magnetic field to flip the electron spin in the 
qubit dot. For a successful injection, the electron with spin up will flip to become spin 
down and then hop into the right dot. The DQD system then takes the state | IlTr)- 
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As discussed above, this state, i.e., the superposition state of |3) and |4), is coupled to 
the single-dot singlet state |5), giving rise to a variation of the occupation probability 
CX55 [see Figure. 0(a)] as well as the QPC current [see Figure. 0(b)]. This is different 
from the case of a constant current in the absence of any successful electron injection 
into the left dot. 

4. Conclusion 

In summary, we have studied the readout of a single electron spin in a DQD system. 
The electron spin is initially confined in the QD serving as a qubit dot. A reference dot 
is coupled to the qubit dot via a tunneling coupling. Also, a QPC acts as a measurement 
device, placed near the reference dot for detecting the variation of the electron number 
in the reference dot. We have considered the two implementations in which either one or 
two electrons occupy the DQD. In the one-electron case, the only electron in the DQD 
is the qubit electron to be measured. An external magnetic field is applied to both dots 
so that the energy level splittings and for spin-up and spin-down electrons are 
different. Gate voltages of the two dots are tuned so that e-f ~ and e± ^> f2 . These 
conditions ensure that only a spin-up electron but not a spin-down electron in the qubit 
dot can tunnel to the reference dot. This gives rise to very different currents through 
the QPC and can be used to readout the electron spin states of the qubit dot. In the 
two-electron case, an additional spin-up electron is always confined in the reference dot. 
This can be easily achieved by properly tuning the gate voltages of the dots. We have 
shown that the electron spin states of the qubit dot can also be readout by applying an 
external magnetic field when considering effects of hyperfine interactions between the 
electron spin and the nuclear spins of the host materials. In the high Zeeman splitting 
limit, the flipping of the electron spin induced by the hyperfine interactions are greatly 
suppressed. In this case, only a spin-down electron in the qubit dot can tunnel to the 
reference dot. This again allows one to distinguish between the electron spin states in 
the qubit dot by measuring the currents through the QPC. Furthermore, we propose an 
approach involving an ESR oscillating magnetic field which can confirm the success of 
an electron injection event into the qubit dot. 
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